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Abstract 



The canonical partition function of a system of rotators (classical X-Y spins) on a 
lattice, coupled by terms decaying as the inverse of their distance to the power a, 
is analytically computed. It is also shown how to compute a rescaling function that 
allows to reduce the model, for any d- dimensional lattice and for any a < d, to the 
mean field (a = 0) model. PACS: 05.20.-y, 05.70.Ce, 05.10.-a 
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1 Introduction 



Let us consider the following classical haniiltonian model of a system of rotators: 

N N 

^=2^^^ +2.^ [l-cos(0,-0,)]=i^ + y. (1) 
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The potential energy V is not thermodynamically stable and the ensemble av- 
eraged energy density U = (^) diverges in the thermodynamic limit (TL) Q. 
If the potential energy term is divided by N , then the energy density becomes 
intensive and it is bounded as N goes to infinity. 

Indeed, dynamics and thermodynamics of the 1/iV rescaled model has been 
extensively investigated in particular, Ruffo and Antoni, who called it the 
hamiltonian mean field X-Y model (HMF), solved it in the canonical ensem- 
ble, and compared the theoretical caloric (T vs U) and magnetization [M vs U) 
curves with those obtained from a microcanonical simulation |3|. 

Here we consider a generalization of model 

ij^iyx^^iyi^f^fcM. (2) 
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The rotators are placed at the sites of a lattice and the interaction between 
rotators i and j decays as the inverse of their distance to the power a. 

A onedimensional version of model (||) has been studied by Anteneodo and 
Tsallis who have numerically measured the largest Lyapounov exponent, as 
a function of N and a. Through a rescaling factor N* — Anteneodo 
and Tsallis showed that their results coincide with those previously obtained 
for the HMF [a = 0) model; this rescaling could then give a well defined TL to 
model @). 

In a recent paper Tamarit and Anteneodo, using a rescaling factor N = 
2" ^ , have shown that the caloric and magnetization curves of model (||) in 
one dimension collapse onto the curves of the HMF model |^ . This universality 
emerges plotting T/N as a function of H/NN and M as a function of H/NN, 
from molecular dynamics simulation of model (^) for different N and a values. 
These authors conjecture that the results they obtained in the onedimensional 
case might be general, valid in any dimension d and for a < d, as suggested also 
in S. 



2 Partition function 

In this work, inspired by 1^ and j5|, we analytically compute the partition 
function of an iV-rescaled model (||) for any d and a < d. In formula (^ij) we 
give the right expression of the rescaling function iV, to obtain universal state 
curves for all lattice models with long range (a < d) interactions. 
Let us now rewrite the rescaled version of Hamiltonian (|^) : 

Ix^ .9 1 >^ l-cos(6l, -61, ) 
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where we have introduced an external magnetic field h — {hx, hy) of modulus 
/i, that makes possible to compute the magnetization. The indexes i, j label the 
sites of a d-dimensional generic lattice; rij is the distance between them, with 
periodic boundary conditions and nearest image convention (the definition of rn 
will be given shortly); a > 0. At each site a classical rotator (X-Y spin) of unit 
momentum of inertia is represented by conjugate canonical coordinates (L^, 0^), 
where the L^'s are angular momenta, and the 9iS € [0, 27r) are the angles of 
rotation on a family of parallel planes, each one defined at each lattice point; x 
and y refer to the components of boldface twodimensional vectors defined over 
these planes. To each lattice site a spin vector 

rxvi ^ {niixTmiy) = {cos 9 i,s\n 9 i) (4) 

is associated, and the total magnetization is given by: 

1 ^ 

M=(Af„M,) = -^m.. (5) 
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Note in the rescaling factor N in front of the potential energy term, now 
written as a free double sum over both indexes. N should be regarded as an un- 
known function of N, a, d and the geometry of the lattice, with the fundamental 
property of making 

an intensive quantity; this guarantees the thermodynamic stability of the po- 
tential. We also note that the sum in (||) is independent of the origin i because 
of periodic conditions. To reproduce the usual HMF it is also necessary that 
N{N,a — 0,d) — N. The constraint i ^ j over the double sum is removed 
defining rf^ = 1/6, a finite number. Since the numerator 1 — cos(0i — 9j) is 
zero for i — j the choice of h is free. The removal of the constraint allows to 

introduce the distance matrix i?^ - = -kr] the diagonalization of such matrix is 

-' "^ij 

the key point to obtain, in the computation of the partition function, known 
integrals in the variables 9i. 

As usual the partition function factorizcs in a kinetic part Zk = ( ^ ) i 

where f3 — l/fc^T, and a potential part Zy- After defining Rij — i^R'ij) 

B — /3h, C — exp ( — 7^ X^ii W ) ' potential part can be written as: 



Zv = C I d^9 exp 
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where fi = x,y. Diagonalizing the symmetric matrix R = (Rij) with the unitary 
matrix U such that R — U'^DU, D — (RiSij), where Ri are the eigenvalues of 
R, we can write the first part of the exponent in (|^) as: 

E {mixRijiriyx + miyRijjnjy) = ^ ("-L-^i + n'^yRi) , (8) 
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where rii^ — Uijnijf^. In order to apply the gaussian transformation: 

/ dze-^+^' a > (9) 



to each term of the sum in the right hand side of (|^), each Ri must be positive. 
The spectrum can be exphcitly computed using a d-dimensional Fourier trans- 
form of matrix R, the eigenvalues being labelled by vectors of the reciprocal 
lattice. These eigenvalues are trivially related to those of matrix R'. A study 
of the spectrum of R' in the limit N ^ oo and for 6 = shows that: when 
a > d each element of the spectrum converges to a finite quantity, the least 
eigenvalues being negative and of order one in modulus; when a < d a, part of 
the spectrum converges to a finite quantity, another part diverges to +oo, at 
most as N. However this last part consists of a fraction of the total number of 
eigenvalues which goes to zero in the limit iV — > oo. The least eigenvalue is still 
negative and of order one in modulus. Then part of the spectrum is negative, 
but it is easily seen that it is shifted by b. Thus calling p the least eigenvalues 
of R' for 6 = and choosing 

b=-p + e e>0 , (10) 

we have that with this 6 the whole spectrum of R' (and therefore that of R) 
becomes positive. Then for each i = 1, ■ ■ ■ , N , ji — x,y we can apply with 
the correspondence a Ri, S ^ rij^, z — > z^^. Performing the integrals over 
variables 9i and using the transformation = '^'^jiUR)ij'^jfi with Jacobian 
2^ det R, we can rewrite the partition function as: 

Z = CZk^ l^^d^^^d^^y (11) 

N 
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where Iq is the zeroth order modified Bessel function. The isolation of the N 
factor in the exponential prepares the object for the use of the saddle point 
method. The quantity in square brackets is intensive. Double sums in the first 
two terms are compensated by R/N — {(3 /2N N)R' and the last sum has 1/iV 
in front of it. The argument of /q is also intensive because involves a term of the 
form Rij — {P/2N) J2j R'lj- If we call f{w) the function in square brackets, 
where w = {^ix, ■ • • , Nx^ ^ly, ■ ' ' i ^ Ny)i then the application of the method 
requires the following three conditions: f{w) admits a stationary point iuq; is 
a simple stationary point, i.e., det Hef\o ^ 0, where Hef\o is the hessian matrix 
of / in wq; the path of integration can be deformed (generally going into C'^^) 
into a path that passes through wq following the steepest descent of f{w) and 
such that f{w) < /{wq) throughout the all path. If the point wq is a maximum 
no deformation is necessary and the method is also called the Laplace method. 
Since, as we show below, wq is indeed a real-valued maximum, we readily obtain 
for the free energy per particle F: 

-I3F = lim — - = lim {-In — - ^ > — 



ij 
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The stationary point wq is given by the vector {'^x, • • • , ^f^^, 5'^ , • • • , ^j^), homo- 
geneous on the lattice sites. Defining ^ = {^x, '^y), its direction is that of B, 
and its modulus ^ is given by the solution of: 



h 



(13) 



with 



A- 



N ^ N 



E- 



(14) 



and where Ii is the first order modified Bessel function. In (|T^ A does not 
depend on i because of the periodic boundary conditions. We note that when 
h = we have infinitely many degenerate solutions, since only the modulus ^' is 
determined. Evaluation of the elements of the hessian matrix at the stationary 
point gives: 

N d^f 
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where we do not give the explicit expression of ff^ty(u'o). As we will see shortly, 
the eigenvalues analysis of the hessian matrix (Hq) shows that the stationary 
point wq is a maximum. Then, Laplace method applies and Eq. ( |T^ ) is valid. 
However, only in the long range case (a < d) the last term in the rightmost side 
of ( p^ is zero; when a > c? its expression does not appear to be manageable. 
We will comment on this point later. Restricting then to a < d, and comput- 
ing the derivative of (|l2|) with respect to the magnetic field we find that the 
magnetization M = (|M|) is given by the solution 5* of (p^). Then the internal 
energy U is given by: 



df3 2/3 2 ^ ' 



(16) 



Equations and (|T|) are the same as those of HMF, as soon as a proper N 
rescaling gives 



N ^ ' N 

3 



= 1 . 



(17) 



Now, from equations (|T^) and (pT[), and calling A„ the eigenvalues of R' , we 
find that, choosing B along one of the coordinate axes, the eigenvalues of the 
hessian matrix at the stationary point are given by: 



(1) _ /? A„ 
2 N 

(2) _ P 

2 N 



* + /i TV 



(18) 



Following our previous analysis we have that: 



e A„ 
N ~ N ~ 



(19) 
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Then we immediately see that Xn^ ^re all positive for any (3 and h] for 

we need to include '^{(3,h) from ([l3|). We have checked numerically that the 



(1) 

n 

are also all positive. From (|18|) we can derive an expression for the determinant 



quantity in round brackets in (llq) is always smaller than 1, and therefore x 



of matrix (15). It is given by: 



1 

N 



In dot ( -yiJe/loj = ^Indeti? 
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N 
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N 



* + /i TV 



}■ 



(20) 



When a < d most of 4?" SO to zero for 



oo, then the sum in ( pO| ) is effectively 
constituted by the terms with the remaining These terms are a fraction of 
N that, as we already pointed out, goes to zero when N oo. If we call N'{N) 
this fraction, then the sum in (^0|) can be bounded from above by ^j-c for 
N oo, where c is a finite number. Therefore the last term in (|l2| ) is zero. 
When a > c? all terms contribute to the sum in (^o|), and we can not give a 
meaningful expression for (|l^). At the end of the calculations we can let e ^ 
in (0). 

Then we have shown that any model with a < d on any lattice is equivalent 
to HMF. From we get an exact expression for N: 



(21) 



We have made a microcanonical simulation of Hamiltonian (^) on a threedimen- 
sional simple cubic lattice in zero magnetic field, using a fourth order simplectic 
algorithm ||^ with time step 0.02, selected to have relative energy fluctuations 
not exceeding 1/10^. We have chosen a fixed N = 343 = 7'^, and have simulated 
various energy densities H/N and various a < 3. In Fig. 1 we show that the 
numerical caloric curves collapse onto the universal HMF curve. The kind of 
results shown in |^ for a onedimensional lattice, where a slightly different N 
has been used. 



3 Conclusions 

Going back to the beginning of our discussion: it is now clear that model (||) 
completely reduces to model (^ for a = 0. In model (|^) the range of the 
interactions is infinite; each rotator interacts with all the others and with the 
same intensity. To get a well defined TL it is sufficient to divide V in 
by N, the total numbers of rotators. It is then possible to compute caloric and 
magnetization curves |^ ; the spatial arrangement of the rotators has no effect on 
them since the intensity of the interaction is the same for each couple of rotators. 
In this work we have shown that, when considering model (^, it is possible to 
take into account the spatial c?-dimensional arrangement of the rotators and the 
decaying of their mutual interaction through a factor N, which is computable 
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Figure 1: The full line gives the canonical theoretical caloric curve (temperature 
T vs energy density U) for long range rotators compared with the microcanon- 
ical simulation of a threedimcnsional simple cubic lattice for three different a 
values: 0.75 (open circles), 1.5 (diamonds) and 2.25 (crosses). Note that in spite 
of the size of the system, still not very large (side with 7 lattice sites) , the results 
already follow very well the theoretical curve. 



for any periodic lattice and any a < d. Dividing by N the potential energy 
in (|^), the model gets a well defined TL and it is possible to compute state 
curves which become those of the HMF model with a proper normalization 
of the constant A in (p7|). The HMF (a = 0) model has revealed peculiar 
equilibrium and nonequilibrium properties , namely: ensemble inequivalence, 
metastability, collective oscillations, anomalous diffusion and interesting chaotic 
properties, both in the ferromagnetic and antiferromagnetic case. On the basis 
of the thermodynamical equivalence here established it would be interesting to 
investigate the a dependence of all these properties. The study of the Lyapounov 
exponents in B is the first in this direction. 
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